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Ce document est un support pédagogique du cours destiné aux étudiants de 2ème Année des licences 
LMD assurées au trois départements: Maintenance en Electromécanique, Maintenance en Instrumentation 
et Hygiène et Sécurité Industrielle aux seins de l’institut de maintenance et de sécurité industrielle. Dans 
ce cours on s'intéresse a un certain nombre de méthodes directes et itératives utilisées pour la 
résolution des systèmes d’équations linéaires, des équations non linéaires, l'approximation des 
fonctions par interpolation polynomiale, valeurs propres et vecteurs propres, la résolution numérique 
des équations différentielle et les formules de quadrature, Trapèze et Simpson pour l’intégration 
numérique. Le tous est regroupé sous le terme générique de “° Méthodes Numériques ‘”. Toutes ces 
techniques de résolutions permettent de résoudre de manière exacte ou de manière approchée un 
problème donné. 


A la fin de ce cours, l’étudient obtient des connaissances solides sur différentes méthodes numériques 
que par la suite sera capable de les implémenter en langage de programmation tel que Matlab ou Scilab. 
De plus ce cours est une base au plusieurs cours de méthodes numérique qui seront enseignés durant le 
cursus de graduation de nos étudiants. 
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Introduction 


Afin de bien maitrisé les méthodes numériques, J’ai pensé de présenter un rappel sur l’algèbre linéaires 

plus précisément sur les matrices (Critère de Cauchy pour la convergence) pour introduire aux 
étudiants les techniques de bases utilisées pour la résolution des systèmes d’équations algébriques. J’ai 
abordé dans ce support quelques méthodes directes et indirectes qui nous permettrons de résoudre les 
systèmes d’équations linéaires tel que (Substitution, Cramer, élimination de Gauss, élimination de 
Gauss-Jordan, inversion de la matrice et la factorisation LU d’une matrice ), des équations non 
linéaires tel que (La dichotomie, La sécante, Newton Raphson), Lagrange comme une des méthodes 
d’interpolation polynomiales et le polynômes caractéristique des matrices pour la résolutions des 
valeurs propres, des vecteurs propres, l’inverse et les puissances d’une matrice. On termine par 
deux méthodes (Euler et Runge -Kutta) pour la résolution numérique des équations différentielle ainsi 
que les formules de quadrature, trapèze et Simpson pour l’intégration numérique. Ces méthodes vont 
êtres implémenté en langages de programmation Matlab, pour cela j’ai présenter une annexe qui résume 
un exemple élémentaire de script qui pourra êtres utile aux étudiants lors des séances de TP et des 
sujets d’examens de synthèses et de rattrapages depuis 2014 au 2018 . 


Le contenu de ce polycopié de cours intitulé «Méthodes Numériques » résumé en sept chapitres suit 
conformément le programme ministériel de l’enseignement supérieur et de la recherche scientifique 
proposé aux étudiants de Ze: année LMD avec un Volume horaire hebdomadaire : 01h30mn de cours et 
01h30mn de TP) ; Coefficient : 02 ; Crédit : 4. 


Les mots clés 


Matrice, Matlab, Gauss, convergence, dichotomie, polynôme, Newton, Jacobi, Gaus Seidel, Lagrange 
Euler, Runge —Kutta, trapèze, Simpson. 
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L'organisation du polycopié 


Le cours des méthodes numériques que j’ai enseigné durant ces 5 dernières années pour les étudiant de 
deuxième année licence et rapporté dans ce polycopie est décomposée en sept chapitres organisés 
comme suit : 

Un rappel sur l’algèbre linéaire est présenté au premier chapitre. J’ai abordé en deuxième chapitre des 
méthodes de résolution des systèmes d’équations linéaires (directe et indirecte). Le troisième chapitre est 
dédié à la résolution des équations non-linéaires. Le quatrième chapitre présentera le polynôme 
caractéristique d’une matrice, les valeurs propres et les vecteurs propres ainsi que l’interpolation 
polynomiale. La résolution numérique des équations différentielle est présentée en cinquième chapitre. Le 
sixième chapitre présentera l’intégration numérique avec les formules de quadrature, trapèze et Simpson. 
Nous terminerons par une annexe en septième chapitre qui résume des sujets d’examens de synthèses et 
de rattrapage et un script en Matlab que j’ai proposé dan notre établissement IMSI. 
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Chapitre 1: Rappel sur l’algèbre linéaire 


Sommaire 


1.1 


1.2 

1.2.1 
1.2.2 
1.2.3 
1.2.4 


1.3 
1.3.1 


1.1 Introduction 


Introduction 


Vecteurs et matrices 

Applications linéaires et matrices 
Opérations sur les matrices 

Quelques Matrices particulières 

Critère de Cauchy pour la convergence 


Arithmétique 
Mesures de l'erreur 


Afin de bien comprendre le cours des méthodes numériques. J’ai présenté dans ce chapitre un rappel 
d’algèbre linéaire, plus précisément sur les matrices 


1.2 Vecteurs et matrices 


1.2.1 Applications linéaires et matrices 


Soient E et F deux espaces vectoriels de dimensions finies respectives n et m. On fixe des bases 
B= (ej ) j=1,...,m et C = (ei ) i=1,...,n respectivement de E et de F. Toute application linéaire 
rE L(E, F )de E—F peut être représentée par un tableau de scalaires de n lignes et m colonnes 


suivant : 
Ou iz -+ 
dzi a2 ` 
A= 
anm dy 


R din 
don A=(a;),Vil<i<n 
On notera Vjl<j<m 
S A nm 


Lorsque on dispose que d'une colonne ou une ligne, on obtient un vecteur colonne Y ou ligne X 
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Yi 
X2 


Y = On notera Y =(y,), Vi, l<i<n 


X ala Xz e x, | On notera X =(x;), Vj,1<j<m 


1.2.2 Opérations sur les matrices 


1.2.2.1 Somme et différence des matrices 


La somme ou la différence entre deux matrices est une opération très simple. Elle est notée + ou — 
tout simplement et on a la définition suivante : 
Soit A et B deux matrices ayant la même taille, alors si l'on a : 


A=(a;)Vil<i<n 
Vi1<j<m 


B = (b;),Vi,1 <i<n 
Vi1<j<m 
Alors 
A+B = (a; +b) 
A-B = (a, -b,) 


1.2.2.2 Déterminant 


A est une matrice carrée, on associe une quantité scalaire appelée déterminant et notée det (A) 
calculée par (1) 


det A= D (-D'a, "Az; t Anj, OÙ jio Joses jn (1.1) 


Représentent une des a (a=n!) permutations possibles des entiers 1,2,..., n 
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Exemple 
ES Le 
Ai iz Miz det A = (-1)°a,,a,,a,, +(-1)'a ite + (—1)°a a343, + (1)? a 24, a33 + 
A= 
RNA (Das + (-1) a, 3433431 


Las dy as | 


Le déterminant d'une matrice diagonale est égal au produit de ses éléments diagonaux 


Exemple 
"non 
det(diag( 10,53) = | o s01= Ta =10*5*3-150 
k=1 
0 0 3 


Une matrice diagonale est inversible si et seulement si son déterminant est non nul, c'est-àdire 

si et seulement si tous ses éléments diagonaux sont non nuls. Dans ce cas, l'inverse d'une 

matrice diagonale est une matrice diagonale où les coefficients diagonaux sont les inverses des 
coefficients diagonaux de la matrice de départ. 


En effet, si : gi 
Lg 0 0 


D = éiser nu sun = | 0 Jm 


noh 


1.2.2.3 Trace 


La trace de la matrice carrée À d’ordre n est (2): 


n 
Tr(A)= > d;; (1.2) 
i=l 
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1.2.2.4 Produit d'un scalaire par une matrice 


Soit une matrice A telle que A= (a) etun scalaire o, on peut multiplier une matrice À par un 


scalaire & 
OO, A), OO, 
Or: Qaa "rr Gi, 
OO. QA, <- Ca, 
CA = (aa) = 


Plusieurs remarques sur cette opération : 


On peut ne pas écrire le point de la multiplication ; ainsi on écrit œ A plutôt que œ A 
Le scalaire s'écrit toujours à la gauche, ainsi on écrit œ A mais surtout pas AO 


De même on écrit A A mais surtout pas A/@ 


Œ'(A+B)= 4 A+CUB 
(4 +B)A= & A+BA 
Œa(B A)=(® DA 
1.2.2.5 Produit d’une matrice par une matrice 


Soit à multiplier deux matrices A et B, pour que cette opération devient réalisable : il faut que A possède 
d’autant de colonnes que B a de lignes. 


Le produit de A [n,m] par B [m,q] est la matrice C[n,q] dont la colonne numéro j est le produit de A par la 
colonne numéro j de B et ce pour chaque numéro de colonne j compris entre 1 et q 


pour Arial Vil£i£n ——  B=(b)Vil<i<m 
Vi1<j<m Vj1< j<q 
m 
> x 
alors C=A*B = Cj; = dx D 
k=1 
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1.2.3 Quelques Matrices particulières 


Il existe plusieurs types de matrices particulières, citons quelque une : 


e Matrice carrée 


Si les espaces E et F sont de même dimension n alors À représente une fonction de E dans E 


r:E—>E la matrice A sera de n lignes et n colonnes, A est dite matrice carrée. 


dii az -+ A 


n 


KEISER 


n 


|an Da +++ Apn 


La ligne de gauche à droite décroissante d'éléments de gu à a, est appelée la diagonale 
principale de la matrice A ainsi que La ligne de droite à gauche décroissante d'éléments de o. 


` 


à d est appelée la diagonale secondaire de la matrice A. 
e Matrices égales 
Deux matrices A(n, m) et B(n, m) sont égales si 
a=b; VI<i<n,1<j<m Unnote A=B 
e Matrice symétrique 
C’est une matrice carrée A(n, n) tq: a. =a. 
e Matrice diagonale 
C’est une matrice carrée tq: AaS AEN AER 
1000 a; =0,Vizj 


Exemple 


050 
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e Matrice scalaire 


C’est une matrice diagonale tq: 4 =A. AER 
ii à 


[5 0 0 


Exemple 


050 


0 05 


e Matrice identité 


C’est une matrice scalaire tq: 4=1 


[100 


Exemple 


e Matrice inverse 
On appelle matrice inverse d’une matrice carrée A, notée AT tq: 


A'A!=AT!.A=-I] 
e Matrice triangulaire inférieure gauche 


La matrice triangulaire inférieure gauche est une matrice carrée L tq: } =0 Vi>i 
Nee 


Gi $ 
L, L, O 1 0 0 
Sc A, Exemple | 410 
[692 

AT A 
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e Matrice triangulaire supérieure droite 


La matrice supérieure droite est une matrice U tq: DU, = 0, Vi> j 
(nu, E um, | 
[107 4 
Uz 
U = Exemple 051 
O Bä 
[0 03 
L U nn J 
e Matrice orthogonale 
C’est une matrice carrée tq: Ar. AA. Ai 
e Matrice à diagonale fortement dominante 
A est dite à DFD si elle est (3) tel que : 
Vi € IX a;| > > a; (1.3) 
| G 
1071 
Exemple A= 
-2 5 -1 

ng 

Wall 

"D >F2+f1 

x B (18 +] 


Donc A est non DFD 


e Matrice augmentée 


Soient la matrice A(n, m) et le vecteur b(n), on appelle matrice augmentée que l’on note: /A, b] la 


matrice obtenue en augmentant À par b, b devient la colonne de A. 


Exemple : 


0 27 15 


A 
—100 1.3 18 


0 27 15 -5] 
LA. b] = 


—10 1.3 18 -6 
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e Matrice transposée 


On appelle matrice transposée B d’une matrice carrée À, toute matrice B tel que Bu e 


On notera PD Ai avec A— de même que le vecteur 


X = | Sera noté X =(x...x,) 


e La matrice nulle 


C’est la matrice non nécessairement carrée dont tous les coefficients sont nuls. On la note On,p ou On p 
dans le cas ou elle possède n lignes et m colonnes. 


Exemple: ©O,, = D OO 
O OO 
1.2.4 Critère de Cauchy pour la convergence 
Soit V une norme matricielle sur M(R). Pour qu’une suite de matrices AU) soit convergente il faut 
et il suffit qu’elle vérifier le critère de Cauchy suivant (4) : 


Quel que soit e> 0, il existe po EN tel que les conditions m e€ N et et p > pọ entraînent 


y(APtM) — AP) <e (1.4) 
1.3 Arithmétique 


1.3.1 Mesures de l'erreur 


L’implémentation des algorithmes conduira inévitablement à des erreurs. Pour mesurer celles-ci, 
nous devons mesurer la distance entre le vecteur représentant la solution exactex = (x1, X2 , X3 
a, Xn ) et le vecteur x = (x1, ko ES EE ) ) représentant la solution approchée. Pour 
cela nous utilisons la longueur de R” (5): 


Ixllo = maxx;| 1<i<n (15) 
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Exemple 
Soit x =(1,—7,2,4) alors 


xl = max|x;| =7 
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Sommaire 
2.1 Introduction 
2.2 Les méthodes directes 
2.2.1 Cramer 
2.2.2 Gauss sans pivotation 


2.2.2.1 Gauss avec pivotation 
2:2:2 Gauss-Jordan 


2.3 Les méthodes itératives 
2.3.1 Point fixe 
2.3.2 Jacobi 


2.3.3 Gauss-Seidel 
2.4 Quelque exemples élémentaires 


2.1 Introduction 


Un système de n équations àn inconnues ` su réelles est appelé système linéaire d’ordre n. 
Ce système se présente sous la forme suivante: 


duu Fatz +... +a, X, =b 


a, Xi FX +... +4, X, =b, 


gab Lëtz kk ët, =b 


n 


Les (a; een sont les coefficients. 
1< j<n 


Les b, 1<i<n sont les termes du second membre du système. 


Les sont les inconnues 


Ce système peut être écrit sous forme matricielle de la manière suivante 
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dir dis sus e an x b, 
Dr ee venue az, || xX b, 
gier Beet Mes rg X3 S A,,X, 
Dire A, e aan UE b, 


A est appelée matrice de transformation à n lignes et n colonnes 
Exemple 


x+3y+4z=50 


3x+5y-4z=2 système linéaire 


4x+7y-2z=31 
xt +2y°+6z=9 


6x° +4y+2z° =1 
1 


3x? +4y2+6z° =0 


système non linéaire 


2.2 Les méthodes directes 


Les méthodes directes sont des algorithmes qui permettent de fournir des solutions exactes si elles 
existent en un nombre fini d’opérations. Ex: méthode de Cramer, Gauss, de Gauss-Jordan, 
Substitution .... 


2.2.1 Cramer 


L’existence et l’unicité de la solution d’un système linéaire dépend de la valeur du déterminant de A. 
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l- Si det (4)#0 La matrice est dite régulière et le système admet une solution unique qui est 
donnée par la méthode de Cramer. 


det A) det(A,) q — det(A,) 
EE E E EE An T 
det(A)?  det(A) det(A) 


X 


Le déterminant det (A,) est obtenu en remplaçant la om" colonne de A par b. 


ai Aiz du, in DN Gin ++. A 


n 


a, daa dan, ai Dz Gau, Azp 


det (A) = 


a 


nl an2 d,3 eean i-1 b, mixe don 


L’inconveignent de cette méthodes qu’elle est couteuse en temps et en espace de stockage. 


2- Si det(A) =0 A est dite matrice singulière et le système est dégénéré. 


Exemple 


Soit le système suivant 


x+3y+4z=50 2 
La matrice A du système est: — — 
3x+5y-4z=2 ge E 
dE GE 
Donc 
1 3 4 
det A=|3 5-4 =-8Z20 Le système est régulier et admet alors la solution unique : 
4 7-2 
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503 4 150 4 1 3 50 
2 5-4 3 2—4 352 
[317 -2 -24 [431 -}_ -40 . pen hag 
RE a EE -8 -8 


La solution est 


x = (3,5,8) 


2.2.2 Gauss sans pivotation 


La méthode d’élimination de gauss est une méthode directe 


d’équation linéaire Ax=b en un nombre fini d'opérations. 


Elle consiste à transformer le système A*x=b 


triangulaire supérieure droite. 


qui permet de résoudre un système 


en un système £’*x=b’, où A’ est une matrice 


La matrice 4’ est obtenue par substitutions successives dans la matrice d’origine selon l'algorithme 


d’élimination de Gauss suivant: 


1) La triangulatisation 


pour k = 1 jusqu'à n —-1 
pivot <— Ou 


si pivot + Q alors 


(“stratégie de pivot *) 


pouri = k +1 jusqu'à n 


dik 


b, <— b, — 


S b, 
pivot 


pour j = k + 1 jusqu'à n 


Fin pour 


Fin pour 
Fin Si 


Fin 
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On appliquant cet algorithme de triangularisation le système devient le suivant 


di ‘°° à; di, X b, 
0 : ; D'où la solution est 
0 0 a; a, x |=| b; b, 
D W; Xn 8 
: Or, LE : : a, 
1 n 
0 0 O0 an X, b, Seet Su 
KÉ j=i+l 


nn 
pouri =n -l1 jusqu'à 1 
somme <— 0 
pour j =i +1 jusqu'à n 
somme <— somme + ax; 
Fin pour 
bi — somme 
Xi 
d;; 


Fin pour 


3) Le déterminant 


La méthode de Gauss permet de calculer le déterminant par (1) 


Dert A1 = Te (2.1) 
i=l 


avec p est le nombre total de permutation de lignes et de colonnes 


det(A) = (—1” II pivot, 


i=l 


Puisque il s’agit de gauss sans pivotation le p = 0, aucune permutation de ligne ou colonne ( le pivot est 


< > de zero) 


det(A) = GU [pivot, = ] [pivot 
i=l i=1 
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2.2.2.1 Gauss avec pivotation 


e Pivotation partielle 


Dans le cas ou un des pivots est nul, on cherche un pivot non nul dans la même colonne. 
Le pivot choisi est le plus grand en valeur absolu, il est appelé pivot maximal par colonne, celui-ci est 


trouvé on permute deux lignes, ce type de pivotation est appelée pivotation partielle (2) 


Le pivot est choisit par la formule suivante : 


ak" =MAXQaK "| i e [k.n] e» 


Exemple : 


x+3y+42 = 50 E 

à On permute les lignes 2 et 3 x+3y+4z = 50 
GE Ox+7y-—2z =31 
re Ox+0y—47=2 


e  Pivotation totale 


Dans le cas de pivotation totale (3), le pivot choisi est le plus grand en valeur absolu dans toute la 


matrice, il est donné par la formule suivante : 


a = MAX 


rm 


(k-1) o 
aj | K<i=<n (2.3) 
K= Ee 


Et on permute les colonnes r et k et les lignes met k 


Exemple 1 : 
4x+771y-—27 =31 11y+4x-2z=31 
x+3y+4z= 50 S y—21=3 y+4x-2z=3 
On permute 3x+5y-4z=2 On permute 5y+3x-4z=2 
3x+5y-4z=2 
La ligne 1 et 3 la colonne let 2 
x+3y+47=50 3y +1x+4z = 50 


4x+77y-—27 =31 


Exemple 2 : 
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Résoudre le système suivant par la méthode de gauss 


2x "FA, 2x = —6 
a. A "At + x = 0 
3x, X 2x, = 8 
pivot (1) S 2x. ET 
ER T S = —0 
p x +3x, + x, = 0 
(2) = (2)-ar/pivot (1) An — X, + x, +21, = 8 
Se AC, Es = 
2% EX. 2% = —6 
OS E E EE Ee 
FR E EE 
- x, SR, +x, =6 
(3) = (3)-azı/pivot (1) 2x "än -2% = 6 
0 +x% Em HE e 3 
0 = 7x, +4x, +2x, = 17 
sand ao ts — a0 


2.2.3 Gauss Jordan 


La méthode de Gauss-Jordan est une méthode très utilisés, elle agit sur les lignes d'une matrice en 
utilisant les mêmes opérations que celles de Gauss. Les opérations autorisées sont: 


1. Échange de deux rangées. 
2. Multipliez n'importe quelle ligne par une constante non nulle. 
3. Remplacer une ligne par la somme de sur la ligne et un multiple constant d'une autre 
ligne. 
On forme la matrice augmentée correspondante au système d'équations linéaires, après 
application l’algorithme de gauss, on obtient la solution X ainsi que la matrice inverse 


A I b>I A” X 


Exemple 


1 3 31l(x, 0 
Soit le système à résoudre : |2 2 Of[x;|=12 
3 3 61]|x; 11 
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On forme tout d’abord la matrice /A, b] : 


K=1 k=1 
3 3 0 F- 3 -3 -Ù 
Normalisation > |2 2 0 2 [| Réduction — |0 -4 -6 2 
3 é 11 ep: 23; TI 
K=2 k=2 
1.3 3 0 LB =ar e 
Normalisation > [0 1 3/2 —-1/2| Réduction — ID 1 3/2 -1/2 
H: =G 3 11 0 0 6 g 
K=3 k=3 
1 0 -3/2 3/2 1 OÙ -712 
Normalisation > [0 1 3/2 -1/2| Réduction — 10 1 0 -5/2 
0 0 1 413 0 0 1 4/3 


x =(7/2 -5/2 4/3) 


l X; =d; : 
La solution est ° intl. dron: 
2.2.3.1 Inverse d’une matrice 
Selon la méthode de Cramer, une matrice À de rang n n’est inversible que si son déterminant À est 
différent de zéro. Dans ce cas, le produit de A par la matrice inverse A” donne la matrice unitaire 7. 
AA AA = Zon Ais A 


En appliquant la méthode de Cramer sur la matrice A, on peut déterminer A”. 


Exemple : 


wa A 1 
Gi Ge en utilisant la méthode de 47 Be 2 1 =j Sch vérifie Aa AA; 
ramer : 2 que : 
Ber, 
3 6 3 


De même L’algorithme de Gauss-Jordan présenté au début de ce cours (méthode du pivot) opère aussi le 
passage de la matrice C=/A,y] à la matrice D=/1,X] où X est la solution du système linéaire A.X=y ; Soit 
X=A"!y. 
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Après les opérations de l’algorithme de Gauss-Jordan, on obtient : D=A".C=A".[A,I]=[LA"] 
Cette méthode, de calcul de l’inverse d’une matrice qui est résumée ci-dessous, permet de calculer A” 
avec un nombre d’opérations nettement inférieur à celui de la méthode de Cramer. 


Transformation LA. =>(LA”) 


Pour Es L&.,# ona 
Sap = E 
Pour Z Zë "Tu 
i=12,...,& 
Pour j = 27, age ës "ëtt 
ZS 
Exemple : 
1-3 
Soit la matrice A=|2 2 0f . Calculer A” par la méthode de Gauss Jordan. 
3 2 
k=1 
UN er 10W H. 3 3 1 0 Q 
Normalisation > |2 2 0 0 1 0f ;Réductin > |0 -4 -6 -2 1 0 
3 2 6 0 0-1 Der E ef 1 
k =2 
l ESS 1 0 0 LH. 32 <li2 té 0 
e 0 1 3/2 1/2 -1/4 0] *. SE A Er A "ert 0 
Dr og, =3 0 1 0: 0 15/2 1/2- 74 1 
k =3 
1 0 -3/2 -1/2 3/4 ; ED. sëtz: SE 1/5 
Te O 1 3/2 1/2 -1/4 o | Réduction |0 1 0 2/5 1/10 -1/5 
0 0 1 1/15 -7/30 2/15] > 0 0 TI =7/30 2115 
= 2 2 1 
"in e us a 


I3. =F 6. 213 
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2.2.3 Factorisation LU 


Il est très facile de résoudre un système triangulaire comme celui du Gauss que celui ordinaire de 


type Ax=b, donc on décompose A sous forme LU avec L matrice triangulaire inférieur et U matrice 
triangulaire supérieure 


A=LU 
liE — = M; k 
i À Ou 
Lv = 
de D Ly 
(2) Ux=y 
1 0 0 
U 0 | 
A M®P=|0 0 1 0 
8 L Link 
| 0 -l 0 1 
On multiplions la matrice A par une matrice M® 
A“ =MPA® A°=A et A”=U bb, - Ep, 
dx 
Ou 
U=M"®.M®..M A = MA SECH 
ge dr E 


A=M "U =LU 
donc L=M 7 


Les matrices élémentaires M" sont inversibles et leurs inverses sont les matrices L" triangulaires 
inférieures telles que : 


fl 
L? =Jsauf Lei i=ln 19 =1-(M® -1) 
sauf l, =l, i=k+ln 
1 0 0 

1 0 " 0 gw Ze me 
U :. Fa 5 à IP =|0 0 1 "0 

M®2=10 0 1 %0 EO us à 
"fruu : O > Lu 0 1 
O —l 0 1 


LD ON oO 


Algorithme LU 
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pourk = 1 jusqu'à n-1 


pivot <- Ou 


si pivot + Q alors 


Lee] 


(“stratégie de pivot *) 


pouri = k +1 jusqu'à n 


lL & 


pivot 


pour j = k +1 jusqu'à n 


a; <a; -Ll kay 


fait 
fait 


sinon érreur 


fait 


23 Les méthodes itératives 


Les méthodes itératives ou indirectes sont des algorithmes qui permettent de fournir des solutions 
approchées_si elles existent en un nombre infini d’opérations. Ex: méthode de Jacobi, Gauss 


Seiedel.... 


2.3.1 La méthode de Jacobi (4) 


Soit D la diagonale de la matrice A, et G le reste : 


A = D+G 


D 


© 
Il 
Se ES 


0 0 0 0O ap Be 
0 0 0 ana 0 ji 
a; O 0 G=|a, Gi 0 
0 0 diin 
0 D an da di 


Dx"™ =b-Gx™" & x" =D” (b-Gx°) 
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Donona Dx"*" =(b-(E+Px*): M = D (2.4) 
N=-E-F={(E+F) 


Théorème de convergence : 


S1 À est une matrice à diagonale strictement dominante, alors la méthode de Jacobi 


converge 
Démonstration | 
"Jacob: 
vii =(D)'(b-(E+F)x*) 
x“! =Cx+d avee C=D'(E+F), des D''b 
i 1 d 
Ki. = |D 7E +P)|_ =max| — äs 
S Kë i=l,n 
FAIT j= 
2%, bk et = —6 
dë +3x, + x, = 0 
Ax =b & 
3x4 -x +n +24 =8 
- x, +2x, + x = 
soit (%,,%,, X,, X, June solution approchée 
Imaginons que AX=b avec parexemple dist(a, À)< E 
—6-4%, +2X, 
K= S 
0-3-7 e E 
si on essay ait: Dub de SO ee 
_ 8—3% +X, 2%; i a; 
E EE 
6+7, -2X 
X, = 2 3 
1 
tant que dert An, bizse  (ie.107°) 
b, — dan 
ee = j=L j#i 
di; 
fin 
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2.3.2 La méthode de Gauss Seidel (5) 
Soit E la triangulaire inférieure et F la supérieure de la matrice A : 


A = D+E+F 


0 0 0 0 O ap i in 
will 0 0 : 0 0 aji : 
E = di d. il 0 0 0 F ES 0 0 0 di; din 
dia 0 0 : 0 0 an,n 
a n ani a, n-l 0 d 0 0 0 
(Op Elon =b- Eet es am 2(D+E) (b-Fx*) 
M =D+E 
(D + E)x"™ =(p- Fx): | (2.5) 
Exemple 
2%; +4x, 2% = —6 
Xi +3x, + x, = 0 
Ax=bS 
Ap, e Aën + Dë +2x, e 8 
= x, +2x, + x, = 
—6—4X, +2X 
x, = 2 3 
2 i-l n 
0-x -X z e > 
x = 1 %4 b, Out: KE 
~ Soit x, = E E 
8-3x +x, — 2x, ' a, 
X3 — ii 
1 
6+x,—2x 
x, = 2 3 
1 
tant que dert An, b) >£ (Gell) 
i—l n 
new old 
b,- ax) — as 
new j=1 j=i+1 
X, = 
di 
fin 
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Lorsqu'une inconnue est utilisée, c’est automatiquement la plus récente valeur qui est calculée. 
Ceci assure une convergence des calculs bien plus rapide que la méthode de JACOBI. 
On arrête les calculs lorsque les valeurs successives de x; sont suffisamment voisines. 
Pour cela, on peut utiliser, 

e soit le critère de Convergence absolue (6) : 


Fil — (E SG 
x; — x; ZZ (2.6) 


e soit le critère de Convergence relative (7): 


(ai) _ IR] 
Ti AE tee e (2.7) 
CET Gg 
xy 


Pour les systèmes où les matrices qui sont de rang élevé, il n’est pas commode de faire le test de 
convergence sur chaque inconnue x;. 

Dans ce cas, on fait le test soit seulement sur certaines inconnues que l'on choisit, soit les quantités 
suivantes : 


1/2 k+l) (E) kal) — Li 2 
D # 2 nxt x; # Le x} 
ikel) L pik) (k+l) ik) Ke en 
>k tj | ou Zb xj | FE 2 men |O" ES de Bos 
= FS 1] Ze 


o 
La convergence du procédé ne dépend pas du choix des valeurs initiales e , mais seulement des valeurs 
des coefficients. 
On montre que la convergence est assurée si on a, pour chaque valeur de i (c’est à dire pour chaque 
ligne), la relation est vérifiée. 


N 
SE Zb, SS 
J-l 
Jei 
Autrement dit, il y a convergence si chaque élément diagonal est supérieur ou égal, en module, à la 
somme des modules des autres éléments de sa ligne (la matrice A soit DFD (8)). 


2.4 Quelque exemples élémentaires 


Exemple 1 : Voici un aperçu d’exécution pour le système suivant par la méthode de Gauss 


4x =2x, 2x = —2 
4x EX 24 = -11 
Ax =b © 
24 +34% e = -1 
4 -2 2 -2 
4-3 -2 -11 
2 3 -1 -1 
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En appliquant l’algorithme de triangularisation de gauss on trouve le système suivant 


Après substitution on trouve la solution suivante 
X=(-1 ,1,2) 
Le déterminantest det(A) = Cl [pivot, =4*-1*-18=72 
i=l 


Exemple 2: Voici un aperçu d’exécution pour un autre système par la méthode de Gauss 


St, +2x, -x =5 

ere x ` deg, e = À 

E 2x +x, +4x, = 7 
5 2 -1 6 
1 6 -3 4 
2 1 4 7 


Après substitution on trouve la solution suivante 


X=(1, 1, 1) 


Le déterminant est det(A) CU [pivot, =121 


i=1 


Exemple 3 : Voici un aperçu d’exécution pour le système suivant par la méthode de Gauss Seidel 
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2%, RE EX = 13 
fe Me. Se VE = 8 

4x, 5 +7x, = 16 
2X1+X2 -X3= 13 X1= (13 -X2+X3)/2 


X1 -5X2 -X3 = 8 > X2=(8 -X1 +X3)/-5 > 
4X1 +0X2 +7X3 =16 X3 =(16 - 4X1)/7 
Après 3 itérations et x(0) = (1,0,0)7 
On trouve la solution suivante x) = (6.0139050, —-0.16744637, —-1.1508029)7 


Exemple 4: Voici un aperçu d’exécution pour le système suivant par la méthode de Gauss Seidel après 
la vérification de la convergence 


Soit le système Áx =b 


2 0 7 x1 16 
-2 5 *| x2] = |23 
6 -2 2 x3 0 


En algèbre linéaire, une matrice carrée à coefficients réels ou complexes est dite à diagonale dominante 
lorsque le module de chaque terme diagonal est supérieur ou égal à la somme des modules des autres 
termes de sa ligne. 


A = (a), Vi,l <i<n 
Vi1<j<m 


n 


a;|> > la 


j=l 
Zei 


on a alors 


Vie D. n|, 


y 


On vérifier la formule pour les trois lignes : 


2] < 0| + |7| pas bon 
|5| > Lü ` bon 
DI < ek pas bon 
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Donc A n’est pas dominante 


Pour avoir A strictement diagonale dominante il faut arranger A, donc si on permute la 1 colonne avec 
la 3 colonne et on change la position des Xi, on obtient le système suivant 


70 2 x3 


16 
0 5 -2| *|x2| =| 23 
2 -2 6 x1 0 


[7] > [0] + [2] bon 
[5] > üb |-2] bon 
16 > Dä ` bon 


Une fois notre système set DFD, on exécute méthode de Gauss Seidel avec 5 itérations et comme 
solution initiale x® = (1,0,0)7 , On trouve la solution suivante 


X3=2, X2=5, Xl=1 


Donc la solution est x) = (1,5,2)7 


Exemple 5: Voici un aperçu d’exécution pour le système suivant par la méthode de Gauss Jordon 


SE, 42%, 63%, =g 
2X, —X, TX; —=8 
3x, 


— X3 = 3 


Etape 1 Former la matrice augmentée correspondant au système d'équations linéaires 


1 2 3 9 


Etape 2: Transformer la matrice augmentée en matrice sous forme d'échelons réduits par des opérations 
de lignes élémentaires. 
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LA 


Étape 3: Résoudre le système linéaire correspondant à la matrice en forme d'échelon de ligne réduite. La 
solution (s) sont également pour le système d'équations linéaires à l'étape 1. 


x = 2, x% = —1, x, =3 


Exercice 


Est ce que le système linéaire A x= b suivant converge et pour quoi ? 


10 2 3 x —1 
a 2 12 1 | ell -{ 2 | 
1 4 10 D 3 


- Donner les matrices D, U et L qui correspondent à la décomposition A =D- U-L dela 
méthode de jacobi 
- Résoudre le système linéaire A x=b avec x(0) = (0,0,0) en 4 itérations. 
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Résolution des systèmes non- 
linéaires (zéros des fonctions) 
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fonctions) 


Sommaire 


3.1 Introduction 


3.2 La dichotomie 
3.3 La sécante 
3.4 Newton Raphson 


3.1 Introduction 


Plusieurs méthodes existent dans la littératures (Point fixe, dichotomie, la sécante, Newton, 
lagrange....) qui cherchent les zeros de fonction = 0, la plupart de ces méthodes sont itératives 


elles calculent des approximations successive is as aen de la racine x de l’équation f(x) = 0, 
à partir d’une valeur initiale x, bien choisie. La seule différence entre eux c’est la vitesse de 


convergence (Figure 3.1). 


- Si f(x) est un polynôme de degré n, on sait que l’équation possède n racines complexes. 

- Si l’équation est transcendante, elle peut avoir un nombre fini, voire nul, ou infini de racines. 
Le problème est alors de trouver la racine dont on sait l’existence et dont, parfois, on connaît 
une valeur approchée. 


D’une manière générale, on suppose que l’équation a été mise sous la forme f(x) = 0 : (ceci est 
toujours possible en définissant par exemple g(x) = x + f(x) puisque lorsque (f(x) = 0, 

g(x) = T). 

x, = f(x) 


x; = f(x) 
x4 =F (x) 


Xy T Zil 
À partir d’une valeur initiale x;, que l’on se donne, on engendre la suite : 


Si la suite des mesures x}, X2, X3,..., x, converge vers une valeur xs, alors : 
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fonctions) 
lim ze bm Zi = Xo ZA) Se : 
#0 Kees „et H ; xo est une racine. 
n= f(x) Kë 
4 A 
~- --------- Ai 3 
1 A K E 
i A A. | n Régime oscillatoire convergent 
1 A k | | 8 7 
w- i l | 
Za ` Sab 
ES wd tee, 
FA | | | | r ai YN 
/ 1 f | 1! i | 
Be ck i | | \/ d 
Lë t 

Y DRE: ge p i | | 

o | w 

: Qui 2. | 
(li E H Ze P OU JS: 2. 7 4 W Lë 


Figure 3.1 : Exemple de solution convergente 


Supposons que l’équation admette une racine x, sur l’intervalle [a,b]. On peut légitimement supposer que 
f(x) prendra des valeurs sur cet intervalle. Si l’on n’ajoute pas d’hypothèses supplémentaires, on ne peut 
être sûr de la convergence. Il est donc impossible de donner une condition nécessaire sans expliciter la 
fonction f. C’est donc une étude de cas. 


3.2 La dichotomie 


La méthode de dichotomie, encore appelée méthode de bissection, est basé sur le théorème de la valeur 
intermédiaire. Elle exploite uniquement le signe de la fonction f aux extrémités des sous-intervalles. 
Elle est décrite comme suit : 

soit, f : [a b] — R, une fonction continue sur l’intervalle [a b]. 


Si f(a)x f(b) < 0 , il existe donc au moins une racine de f(x) appartenant à 
l’intervalle [a b]. 
On prend c =(a + b)/2 tel que : 
1.Si f(c) = 0 — c est la racine de f(x). 
2.Sinon, nous testons le signe de f(a)x f(c) (et de f(c)x f(b)). 
3.Si f(a)x f(c) < 0 — la racine se trouve dans l’intervalle [a c] qui est 
la moitié de [a b]. 
4. Si f(c)x f(b) < 0 —— la racine se trouve dans l’intervalle [c b] qui 
est la moitié de [a b]. 


Ce processus de division, par deux, de l’intervalle (à chaque itération on divise 
l’intervalle par deux) de la fonction est réitéré jusqu’à la convergence (Figure 3.2). 
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Algorithme 


Mé thode de la dichotomie 


a+b 


f(@)fb)<0 > c= 


Tantquela S bé 


faire 
silf (o| <e 
alors on a trouvé la solution: c 
© 
sinon si f(a)f(c) > 0 T 
alors a 4+ c 


sinon si f(b)f(c)>0 
alors b <+ c 


FinTantque 


Figure 3.2 : Méthode de la dichotomie 


3.3 La sécante 


L’idée de cette méthode est très simple : pour une fonction f continue sur un intervalle [a,b], et vérifiant 
f (a) < 0, f (b) > 0, on trace le segment [AB] où A = (a, f (a)) et B = (b, f (b)). Si le segment reste au- 
dessus du graphe de f alors la fonction s’annule sur l’intervalle Ia" ,b] où (a° ,0) est le point d’intersection 
de la droite (AB) avec l’axe des abscisses. La droite (AB) s'appelle la sécante. On recommence en 
partant maintenant de l’intervalle Ta" ,b] pour obtenir une valeur a" ( Figure 3.3) 
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Figure 3.3 : Méthode de la sécante 


Soit f : [a,b] — R une fonction continue, strictement croissante et convexe telle que f (a) < 0, f (b) > 0. 
Alors la suite définie par (1) 


b— an 


-Foran Cw Se 


On = à et unt An 


est croissante et converge vers la solution l de ( f (x) = 0). 
3.4 Newton Raphson 


Lorsque la fonction est différentiable, on peut établir une méthode en exploitant les valeurs de la fonction 
f et de ses dérivées (2). 


Soit x; une valeur approchée de la racine s inconnue. Posons ` x2=x,+h, et cherchons 
l’accroissement qu’il faut donner à zx, de façon à ce que : 


Zei fu +h)=0 
Après développement en série de Taylor à l’ordre 2, on obtient : 
1 A" H 
JG th) =f(Gm)tA f(x) Weed (x + &h) = 0 


ou approximativement ` J) +h (0) = 0 


f(x) 
FX) 


c’est à dire : A=- 


et plus généralement, la solution : **+1 ` 
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(3.2) 


Interprétation géométrique : 


La valeur x; est l’abscisse du point d’intersection avec l’axe ox de la tangente au graphe y=f(x) 
en xy. 


Sens de l'approximation : 


Si l’on avait fait aucune approximation dans l’écriture de Fa +h)=0 


la racine s, l'expression suivante : 


, on aurait obtenu, pour 


Zei KH S'Gt éh) 
“HAL D TG 


ZS" A 


donc : 
o JG) Ai J" + eh) 
kt d WE e 
EE dE 
Si f'(x + 4h) 
Bn, BI ee 
2 T'a) 
b-a 
=b b 
d EE 
c=b- f0) 
f'(b) 
Lu —X— JG) 
TG) 


Les critères d’arrêts peuvent être 
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Es EA Lë 


D xx, 
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Les polynômes : 


Interpolations et Caractéristiques 
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Sommaire 
4.1 Introduction 
4.2 Interpolation et approximation 
4.1.2 Polynôme caractéristique 


4.1.2.1 Trace, déterminants 
4.1.2.2 Théorème de Caylay-Hamilton, 
4.1.2.3 Valeurs propres et vecteurs propres 


4.1 Introduction 


L’illustration la plus triviale dans les problèmes numériques en estimation de données en des 
points intermédiaires non évalués est l’interpolation linéaire où les points de mesure sont reliés par 
des segments de droites. 


On substitue très souvent une fonction f(x) connue en un nombre fini de points x par une fonction 
P(x) simple: c’est l’approximation. En termes mathématiques, l’approximation consiste à 
minimiser la distance qui sépare les fonctions f(x) et P(x). 


4.2 Interpolation et approximation 


Soit f une fonction réelle définie sur un intervalle [a, b]. De par la technologie, l’’etude numérique 
de f ne peut se faire sur un spectre continu de valeurs x € [a, b], mais uniquement sur un certain 
nombre (n + 1) fini de valeurs discrètes (même si (n + 1) peut être aussi grand que l’on veut). C’est 
la discretisation de f 


l’interpolation ` on va chercher une fonction A. passant par y en x pour 0 = i = n et on va ensuite 
travailler avec d 


l’approximation ` on va chercher une fonction 8 qui va d'ecrire globalement le « mouvement » de f 
. Plus exactement, on va chercher ‘a minimiser la distance entre 0 et les (yi) 
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Illustration : 


Les points bleus (e) correspondent aux valeurs connues de f. 


Figure 4.2: Interpolation 


Étant donné un ensemble de doublets numériques (résultats expérimentaux, par exemple), le problème à 
résoudre consiste à trouver un modèle mathématique (polynomial, trigonométrique, exponentiel, etc.), et 
ses paramètres significatifs (c'est à dire ses coefficients), afin de réduire (on parle de régression) toute une 
information en une expression mathématique utilisable, c’est à dire calculable, intégrable, dérivable, etc. 
Lorsque les doublets sont considérés comme ‘sûrs’, au sens expérimental du mot, on tentera une 
interpolation qui restituera toutes les valeurs numériques des doublets là où ils se trouvent. Lorsque les 
doublets sont entachés d’incertitudes sur leurs déterminations, en particulier s’ils sont très nombreux, on 
tentera une approximation qui restituera ‘au mieux’ l’mformation contenue dans les doublets. On raisonne 
sur une fonction numérique ‘f à une seule variable réelle x, connue pour N valeurs. Soit n, le nombre de 
paramètres du modèle mathématique à déterminer. 
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4.2.1 But de l’interpolation 


Étant donnée une fonction f définie sur un intervalle fermé /a,b], pour n+1 valeurs de sa variable x; 
(i=0,1,2, …, n), mais inconnue ou d’expression très complexe en dehors de ces valeurs ; il s’agit de 


calculer une fonction p(z) numériquement plus commode à manier et qui coïncide avec la fonction pour 
les valeurs connues de f. la recherche peut s’étendre au cas où la fonction f n’est même pas connue pour 
ses n+1 valeurs de définition, mais où sont connues des valeurs de ses dérivées, par exemple. La fonction 


; Lu , LCA a E 
P est pratiquement une somme finie de n+1 fonctions de base “i (3), linéairement indépendantes (1): 


j= 


s (2) = D gail E 


Zell 
DICH EE 


base devront pouvoir se prêter à tous les traitements numériques courants. 


le but de l’interpolation est de trouver les a; pour que : et où les fonctions de 


Pratiquement, on pourra choisir la suite des monômes ` 


(x) =] 
u,(xX) =x 
u (x) = x 
t(x) = x 
u (x) = x" 


puisqu’on en vertu du théorème de WEIERSTRAUSS, toute fonction continue peut être approchée 
uniformément par un polynôme. 

Il pourrai être opportun pour simplifier les calculs, de choisir des fonctions de base qui en plus de 
l’indépendance linéaire, ont des propriétés d’orthogonalité, suivant la définition d’un produit scalaire. 
Enfin, pour une fonction périodique, la base formée par les fonctions sinus et cosinus paraît tout à fait 
adaptée. 
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4.2.2 Interpolation polynomiale de Lagrange 


La méthode est ancienne, mais elle est parfaitement adaptée aux traitements informatiques. 
Soit f une fonction définie sur un fermé /a,b] de Di (a<b), pour n+1 valeurs distinctes de sa variable, cet 


ensemble étant désigné sous le nom partition (x; 0<i£x ). 


On cherche s’il existe un polynôme P(x), tel que Peg) = JG) , pour i variant de 0 à n+1 


bad 
= J M 
R(x) KER où les ao, ai, än Am © Di 
Zell 


On a donc : 
mn E 58 AE o a a TE Zu 
bon A aa o ` GE A IEN JD 
E ET A mp E e, E EN "Lë Ge } 
2 3 j 
KE E ff EE S E ët 1 A Sé Se ) 


Dans ce cas, trois éventualités peuvent se présenter : 

e m>n ®© système impossible [| aucune solution, 

e m=n = système de Cramer [1 une solution unique, 

e m<n système surdéterminé O une solution est à chercher par la méthode des moindres carrés. 
Ici, c’est le cas où m=n qui nous intéresse. La solution existe et est unique, car le déterminant des 


Posons : 


"'S 


# 


J 
coefficients ( ZS ) est un déterminant de Van der Monde, non nul et qui vaut : 
A= (xT) 
Otis jón 
Recherche de la solution : 


Considérons les polynômes (de Lagrange) suivants (2) : 


x 


Lœ =| [| (42) 


7 An" 
; ù i=0 àn. 


Le produit effectué sur les indices j tels que TE tin 
DER al "H. ais) 
est clair que et que : 
i a L(x) = 
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L (x;) = d 


Donc, ces  polynômes sont de degré n, et sont tels que 
Ils sont linéairement indépendants, puisque si : 


D 2,1,(x) = 0 


PACE 


yay $ 
alors pour : k aveck=0àn,ona: +? 


Les n+1 polynômes de Lagrange forment une base de l’espace vectoriel des polynômes de degré au plus 
égal à n. Sur cette base, on a : 


P(x) = D f(x) (x) 
i D 
Car, on a bien : R(x) p f(x) pour i=0 à n. 


Autres propriétés : 


=> x Lx) 


Al 
et si l’on pose : 
[I&-x9-[1& 
Jo 
on a alors : 
E E 
(x x)| | (x) 


Cas particuliers des abscisses équivalentes : 


Soit h un réel positif, appelé pas de la partition. 


X; = Xo tik Ze fb tSb 


On pose : `? où s est un réel quelconque. 


Alors l’expression de L (z) donnée dans le polynôme de Lagrange devient (3): 


TXT 2%; les] 
L(x)= Ms, SE Ge Mei (43) 


est un polynôme de degré n, à coefficients purement numériques ( ‘) 


Le tableau des coefficients `Y forme une matrice, de rang n+1, appelée matrice régulière de Lagrange. 
Ces matrices peuvent être tabulées. 


Exemple : Établir la matrice régulière correspondant au cas de l’interpolation quadratique (matrice 
régulière de Lagrange) 
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Solution 


Soit les données de la fonction f suivante : 


x 0 1 2 


zial 1 2 d 


a- SEN 2 
où: 
Recherchons tout d’abord le polynôme d’interpolation de Lagrange de la fonction f . Ce polynôme s’écrit 
sous la forme : 


Aoa D 
SE E EE 
i 0 0 1 1 2 2 
j 1 2 0 2 0 1 


On obtient alors 3 polynômes L (x) i L(x) ; L(x) et L(x) ; 
Laj- hi -3.x+2) 
Lx) = -x7 +3x+0 


L(x) gif - x+0) 


D'où : 
à 2 
Pa) =S LSx) = 2x -5.x+1= f(x) 
bes 
Détermination de la matrice régulière correspondante : 
Ona: 
XTX z SI 
L(x) = i — |= Mis) 
JOUA AS) A 
Ja Ja 
H 
M(s)= zm, Ca | 
Avec j- EIER 


On obtient alors : 


M, (s) = gf -3.5+2) 
Mis) =-5? +35 +0 
Kuel: ZK -s+0) 


D'où la matrice régulière de Lagrange correspondant à l’interpolation quadratique est la suivante : 
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L =3/2 Ms 
Mani 2 ol 
0 -1/2 1/2 


On peut monter également les autres cas de matrices régulières d'ordre inférieur ou supérieur 
(interpolation linéaire, cubique, etc.) : 


1 -1 
M ygi = k i | O interpolation linéaire 


1 =11/6 1 -1/6 
0 3 =3f2 Ma 
Ms = 0 -3/2 > -1/2 O interpolation cubique 
D: 1/3 -1/2 1/6 
L =25/12 35/24 =5/12 1/24 
0 4 —13/3 3/2 -1/6 
Maga = 10 713/12 19/4 S2 1/4 
0 3/2 =713 TI. 1/6 
0 —-1/6 11/24 -1/4 1/24 


Suivant les auteurs et/ou suivant la conduite des calculs, les matrices régulières de Lagrange peuvent être 
ces matrices ou leurs transposées. 


4.2.3 Erreur de l’interpolation de Lagrange 


Pour ZS (xa, äer Ka), il s’agit d’estimer 4x) = d CIE 


l STEE? ; X;) = P (x on 
Si on ne possède d’aucun renseignement sur f, autre que A i) x ( d , On ne peut rien dire sur £ 


dx) =Û0, ap - à 
(autre que éi ). Si l’on ajoute des hypothèses sur f, elles se répercuteront sur €. 


Exercice : Déterminer le polynôme d’interpolation de Lagrange satisfaisant au tableau de valeurs 
suivant : 
x| 0 2 3 5 


y|-1 2 9 87 
Comme il y a quatre points d'appui, le polynôme sera de degré 3. 
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e (x-x)(x-x)(x-x) _(*-2)(x-3)(x*-5)__1 ee 
3 (xo = )(Xo zs )(X0 733) (-2)(-3)(-5) el SE Se 

` Iezllszsstecal a35) Au, 
WNRECAIEEEAIEEES (2)(-1)(-3) d j We: 
pP ee (x° -7x° +10x) 

(% = Xo (X2 -31 )( %2 -%) (3)(1)(-2) 

aeaa) DE- a e 

j (3 =x )(% -4 )(% -%) (5)(3)(2) 50! i Re 


| 53 4 „2 253 
Donc P(x)= y4 (x)+ n4 (x)+ y4 (x)+ 334 (x)= =x -7x + | 


30 30. 


4.3  Polynôme caractéristique d’une matrice 


Le polynôme caractéristique d’une matrice carrée A noté P(A) estle det( A—2AD (c’est un 
polynôme en À) 


| À SE a b a —À b | 8 SET 
Exemple : Le polynôme caractéristique de ( 4) est | S REI ka (A—A) (d-À)-bc 
=2? — (a + d)A + ad — bc 
4.3.1 La trace et le déterminant 


Définition. On appelle la trace de A la somme des éléments sur la diagonale. 


Exemple Soit A une matrice carrée d’ordre 2 


Rappel 
Les valeurs propres d’une matrice carrée sont les racines de son polynôme caractéristique. 


Théorème. La trace de A est égale à la somme des valeurs propres de A et le déterminant de A est le 
produit des valeurs propres de A. 


Démonstration Pour une matrice d’ordre 2 E J le déterminant est 
det(A)= ad — bc, tr (A)=a+d 


det (A-I) = 2? — (a + d)A + ad — bc (1) 
= (à — A)Q À) = 22 —lÀ = A2 + A142 
= 22 = (4 + A )À+ A2 (2) 


De (1)et(2) ona A+ = a+d=tr(A) 
A2 = ad — bc = det(A) 
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Et d’une manière similaire on démontre le cas général 


4.3.2 Théorème de Caylay-Hamilton, 


Théorème : Pour le polynôme caractéristique d’une matrice A, si l’on substitue À par la matrice A, on 
obtient une expression matricielle qui est la matrice des zéros. 


Exemple 
Soit une matrice d’ordre 2 E Si alors 
3 4 


P(A)estle det (A-A) = 12 —tr(A)A+ detif4) = 242 — 54-2 
Selon le thèoreme de caylay Hamilton A? — 54 — 2] doit etre la matrice des zéro, c’est adire 


A-5A-21=0 
Ce théorème nous aide a calculer la matrice inverse ainsi que les puissance des matrice 


La matrice inverse A`! 


A—5A—21=0 
A —5A=2Iet 
A(A-—5)=21et 
AS(A—5)=1 donc 


AT => (4-51) 
Les puissances _ A” 
4 = A2.A=(54+21).A=5 À? +2.A.I = 5(54 + 21) + 2.A.I 


25A + 101 + 2AI = 27A + 101 donc 
A3 = 27A + 101 


De la même manière, on trouve 


Ai = 145A + 54] 


4.3.3 Valeurs propres et vecteurs propres 


Définition 
On dit que å € k est une valeur propre de A s’il existe un vecteur y #0 tel que 


AY= AN 


Dans ce cas on dit que Y est un vecteur propre associé à la valeur propre 1 et que (4, y ) est un 
couple propre de A 


Exemple 


Calculer les valeurs propres et les vecteurs propres de la matrice suivante A= G ` 
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Le polynôme caractéristique de A est P(A)= det(A—Al)= 2? — tr(A)A + defi 
22—-71++6=(1-1)(1—6) 


On obtient donc 2 valeurs propres réelles 1=1 et 1=6 


- On détermine les vecteurs propres associés à À = 1 


Soit le vecteur propre associé vc, 


AY= AN 


E "RI At Sg ES Ke E E EE e KE +2y2= 0 
2 2/7 \Y2 y2 2ÿ1 + 2y, = ds Zu +2Y2 = y2 2y +y: = 0 


a (5) 


=> y, = Än => Y= l 
ae R 


- On détermine les vecteurs propres associés à A= 6 


à Sg yı 
Soit le vecteur propre associé y= 


y2 
AY= AN 

E "FËZk e) E E F2 = rx ES +272 = Din _ E +272 = 0 
2 2/ W2 y2 Zu + 2y, = ds Zu + Aus = Du: 2y1 — Ans = 0 

2 

=> y2=2h =>> y= k Eé 

ae R 
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Résolution numérique des équations 
différentielles 
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Sommaire 


5.1 Introduction 

5.2 Méthodes explicites 

5.2.1 Méthode d’Euler 

5.2.2 Méthode de Runge-Kutta 


5.3 Pas variable 


5.1 Introduction 


Les méthodes de résolution numériques des équations différentielles sont basées sur un 
développement en série de Taylor. 


Pour les équations différentielles ordinaires le lien entre la dérivée et la valeur de la fonction est 
donné par la formule suivante (1): 


ZG FLUG (5.1) 


Avec la forme génirique suivante 


dX(t) 
di = (X().r) — _treprésente le temps 
— _f variable en fonction du temps 
X:R——R" — Parfois Y au lieu de X, parfois x au lieu 
f:R'XxR R’ der... 
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5.2. Méthodes explicites 
5.2.1 Méthode d’Euler 


C’est la plus simple et la plus intuitive des méthodes explicites : c’est une approximation linéaire 
par morceaux de la trajectoire sur l’intervalle [a b], on choisit n+1 points tọ, t1, t1, t2 .... tn, avec 
to =a, tı =b etle pas h = (b — a)/2 


Xo=X(to), on avance d’un pas h avec une taille petite ce qui permet de suivre la courbe de plus 
près (2): 


t =t +h 


(5.2) 
X,=X,+hf(X,,t) 


Dans la méthode d’Euler 


- Lorsqu'on suit la tangente, on s’éloigne de la courbe 


Exemple : f(X,f) = C J 
x 


- Formation d'un cercle lorsque la solution est exacte : 


rcos(t+k) 
Exemple : X(t) = | 
rsint +k) 
- La solution exacte est une exponentielle 
Exemple ` x(t) Ae " 


- _ En fonction de la taille 


x = X(1—-kh) 
h<1/k ok 
Exemple : Gët 
h > 1/k oscillations + / - 
h > 2/k divergence 


- Plus kest grand, plus A doit être petit 
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5.2.2 Méthode de Runge-Kutta 


C’est le même principe que la méthode d’Euler à des ordres plus élevés "Ordre 4 ‘”, cette méthode 
permettre d’obtenir une plus grande précision (3). 


f= TUäaial 
fiz FX, +i fotot?) 
h= e, E 
= FX, +hf,,1, +h) 
Zi Lis X,+4(f+2f+2$, + f,) (5.3) 


5.3 Pas variable 

Com ment choisir le pas h ? 
— Trop large : erreurs, instabilité, divergence... 
— Trop petit : on n’avance pas, long temps de calcul 

e On veut un pas idéal : 
— Aussi grand que possible sans trop d’erreur 
— Lié aux raideurs des équations 
— Le pas idéal peut varier au cours du temps 
L’estimation de l’erreur se fait par deux méthodes : 

— Euler avec un pas h 
— Euler avec deux pas A/2 
— Erreur estimée = différence des deux valeurs : 


Erreur=|X,-X}| 
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6.1 Introduction 


: ER b S Se 
Le but de ce chapitre est d’approcher numériquement f, f(x)dx , très souvent le calcul explicite de 


l'intégrale d’une fonction f continue sur Io, b] dans R, définie par  1(f) = Ip f(x)dx est 
impossible à atteindre. Par conséquent, on fait appel à des méthodes numériques, afin de calculer 
une approximation de (f). Dans ces méthodes numériques, la fonction f, est remplacée par une 
somme finie constituée de n sous-intervalles selon (1): 


b b — 2 
| nie = 2. > feu) 
=1 


Dans ce type d’évaluation, on calcul forcément une approximation (passage d’une intégrale à une 
somme) de la vraie valeur. La méthode d’intégration mise en œuvre est dite d’ordre p si l’erreur 
d'intégration (1): 


Prett) = [fe tie - E2 Rafe) (6. 


est nulle quand f est un polynôme de degré inférieur ou égal à p et non nulle pour au moins un 
polynôme de degré supérieur ou égal à p + 1, soit f € C ([a, b]). Nous allons étudier quelques 
méthodes d’intégration dans lesquelles la fonction à intégrer est substituée par un polynôme 
d’interpolation sur chaque intervalle élémentaire 


6.2 Formule de quadrature 


Soit f une fonction continue sur Io, b] dans R et ont veut approcher numériquement 1(f) = 
b à x 

L f(x)dx , pour passer de l’intervalle sur [x;, x;,1] en [-1, 1] un changement de variable àt est 

nécessaire (2). 
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D'où 


L rd = D D f dx = DAS, rte + Rat (6.2) 


y=f(x) 

N 
h=(b-aXN 
xi=a + ih 


i=0,1,2,....,N 


a b 
x0 x1 xi bi xN 
ZA 

h 


Figure 6.1: Quadrature 


Donc on se retrouve avec une fonction g continue sur [-1 1] dans R et ont veut approcher 
ne 1 
numériquement 1(g) = f; g(t)dt 


Soit M en entier positif M = (1, 2, 3...) et des points d’intégrations comprissent dans l’intervalle 
[-1 1] tel que —1<t <t <t3..tm <1, 


L'idée c’est d'évaluer la fonction en ces points avec des poids d'intégration pour chaque 
point wu W …Wy 


La formule de quadrature pour rapprocher numériquement I(g) = Le g(t)dt est (3) 


J(g) = w1g( t1) + w2g(t2) + +... wug(tm) (6.2) 
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y=g(t) 


1 tt 1 


Figure 6.2 : Points et poids d'intégrations 


Exemple 

On prend un seul point dans [-1 1] 

M=1 ,t4 =0 J(g9)= 2*g(0) 

La problématique est comment choisir les point et les poids d’intégrations 

6.3 Formule du trapèze 

On subdivise l'intervalle [a,b] en sous- intervalles Jl , i = 1,2,..., N; Xo = a; Xn = b} sur 
lesquels la fonction f est remplacée par le segment de droite qui joint les points (x;1, f(x;1)) 
et (x, fx). 

Sur [a,b], on remplace f par une fonction d’interpolation linéaire par morceaux. D’un point de 


vue géométrique, on assimile laire comprise entre le graphe de f et l’axe des x à la somme des 
aires de n trapèzes. 


Méthode des trapèzes 


Figure 6.3: Formule de Trapèze 
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Considérons que la division en sous-intervalle est uniforme et posons : 


X=at+ih où n et f(x) =f; ; i=0, 1,2, ..., n. 


[” Tt 


Sur l’intervalle [xi , xı] Paire est remplacée par h( fi + fi:1 )/2 , aire du trapèze 


correspondant. 


En additionnant les aires des n trapèzes, on obtient la formule des trapèzes (3): 


NOTÉE EA (6.3) 


On peut montrer que l’erreur commise est proportionnelle à h? (si la fonction f est deux fois 
continûment dérivable). 


On dit que la méthode des trapèzes est d'ordre 2. 
La formule est « exacte » pour les fonctions f de degré < 1. 


6.4 Formule de Simpson 


Développons à présent une méthode d'ordre 4 qui équivaut à remplacer la fonction à intégrer par 
des paraboles définies sur des sous-intervalles comprenant trois abscisses d’intégration 
successives. 


Approximation parabolique par morceaux 


Méthode de Simpson 


En noir ‘fonction à intégrer ` ; 


En rouge : approximation parabolique par morceaux 


Figure 6.3: Formule de Simpson 
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On suppose que l’intervalle [a,b] est partagé en n sous-intervalles égaux : 

Ir, x] tels que x;=a+ih, avec h= (Dain, 

On groupe les points par trois, n doit donc être pair 

a = Xo, X1, X2 | X2 X3, Xa | -e |Xn-2 Xn-1 Xn = b. 

Et on remplace, sur chaque intervalle [x;1, x;,,], la fonction f par une parabole. 


Pour l’intervalle [xo, x2], la courbe représentée par f(x) est approchée par la parabole d’équation : 


p(x) = f, (xx )(x— x,) (x- x,)(x — x, ) l (xx LE E 
Lt, z x) — Käl (x z x,)Cx, E x,) (x _ x) E xX) 
o p aaa) p Anaa) (x—x,)(x — x,) 
f 2h? i -h t 2h? 


x2 D h 
LU nee f? pade = SL +4f + F1] 


et l’intégrale est alors approchée par : 
En répétant ce procédé pour les n (pair !) sous-intervalles, on a finalement la 


formule de Simpson (4) : 


oe ELA +45 +2. +4f +..+2f,+4f.+f] (6.4) 


Cette formule est exacte pour les polynômes f(x) de degré < 3 
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Examen de Rattrapage 
Méthode Numériques 
2014 


Exercice 1(5pts) 


Résoudre par la méthodede Gauss, en donnant l'expression de toutesles matrices et seconds 
membres intermédiaires, le système linéaire S1 


2x-y+47 =8 

4X — y +5z +w = 16 

ÉD ee os yes 
3y — 9z + 4w = 3 


(x, y,z,w ) € R* 


- Calculer le déterminant. 


Exercice 2(6pts) 


Donner les matrices D, U et L qui correspondent à la décomposition A = D — U — L de la 
méthode deJacobi dans le cas suivant : 


9 -2 -2 
A=|2 2 8 
3 —9 2 


- Résoudre le système linéaire AX=b par la méthode de Jacobi, b = (11,13,15) 
x® = (0,0,0) 
En 4 itérations 


- Etudier la convergence de la méthode de Jacobi dans le cas suivant: 
2 —1 -1 
A=ļ| 2 2 2 
-1 -1 2 
Exercice 3(5pts) 


-Appliquer la méthode de dichotomie à la fonctionf: x — e fm — xsur l'intervalle IO. 1], avec 
une précision de 10 SS 
Exercice 4(4pts) 


- Donner les instructions Matlab permettant de calculer l’expression suivante : 


5 
1x] + 1 (in(V2x + T)) K 
tg(x?) + (cos(exp x?) + 11) 


- L'écriture A=[2:5 ; 12 :-3 :3] permet d’afficher sur écran: 
A=? 
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Examen Final 
Méthode Numériques 
2015 
Exercice 1 (7pts) 


Donner les matrices D, U et L qui correspondent à la décomposition A = D — U — L de la 
méthode de Jacobi dans le cas suivant : 


9 -2 -2 
A=|2 2 8 
3 —9 2 


- Résoudre le système linéaire AX=b par la méthode de Jacobi, b = (11,13,15)° 
x(0) = (0,0,0)t 
En d'itérations 


- Arranger le système B suivant de manière à avoir une matrice strictement diagonalement 
dominante 
2 0 7 x 16 
B -2 5 Ui D =| 23 
6 —2 2 Z 0 


- Donner les instructions Matlab permettant de calculer les expressions suivantes : 


VRTA(m(RT)) 


arctg (x2)+(sin (exp x?)+1) 


Exercice 2_(5Spts) 


1) 


2) A=[12345;23451;34512;45123;51234 J; 

B = A(1:4,2:3) ? 

C=A :,[1 3] ) ? 
3) Donner des instructions Matlab les plus simples possibles pour produire la matrice A de genre 3 x 3 
ayant la forme suivante : 


Exercice 3 (6pts) 

-Ecrire un script Matlab pour résoudre le système d’équation linéaire par la méthode de gauss 
sans pivotions avec le calcul du déterminant, donner les résultats détaillé 

pour ce système 


5 6 
1 6 3 = 
2 1 4 7 


Exercice 4 (2pts) 

- Citer la différence entre les méthodes directes et indirectes pour résoudre le système 
d’équation linéaire AX=B ? 

- Citez les fenêtres de l’espace de travail du Matlab et sur quoi est basé sans calcul ? 
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Examen de Rattrapage 
Méthode Numériques 
2015 


Exercice 1 (8pts) 


-Résoudre par la méthode de Gauss le système linéaire A x= b suivant et déduire le déterminant 


5 2 —1 x 6 
A=|1 6 -3 ll b=|4 
2 1 4 z 7 


- E = R3 est rapporté à sa base canonique e = {el , e2 , e3 }. f est l’endomorphisme de matrice M dans 
cette base : 
3 1 1 
M = (= 1 al 
0 0 1 
- Déterminer son polynôme caractéristique. 


- Calculer ses valeurs propres. 


- Donner une base des sous-espaces propres. 


Exercice 2 (6pts) 


Décrire l’algorithme de la dichotomie et utiliser le pour calculer le zéro de la fonction 


Tt = x$-4x-8,95 dans l'intervalle [2;3] avec une précision de 0,01 


Exercice 3 Questions de cours (6pts) 
Répondre aux questions suivantes 
- Le choix d’une des méthodes de résolution des systèmes non linéaires est conditionné 
par quels critères ? 
- Donner des instructions Matlab les plus simples possibles pour produire la matrice A 


de genre 3 x 3 ayant la forme suivante : 


- Quand une matrice est- inversible ? 
- Quelle est la condition pour inverser une matrice diagonale ? 
- Pour un système d'ordre 3 comment le résoudre par la méthode de cramer ? 


3 1 1 
M = |-1 1 -1 
0 0 1 


- Déterminer son polynôme caractéristique. 


- Calculer ses valeurs propres. 
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Examen Final 
Méthode Numériques 


Mai 2016 
Exercice 1 (6pts) 
1 Lo Lé La 3/2 
l4 1 14 La 3/2 
Soit la matrice A= | 4 14 1 Lä et le vecteur b = |O 
Lä 4 4 1 94 


- Etudier la convergence du système. 


- Résoudre par la méthode de Gauss en détail le système linéaire suivant et déduire le 


déterminant. 
Exercice 2 (7pts) 


-Pour a € R donné, on considère la matrice M 


2 a 0 
M= |a 2 a 
0 a 2 


-Ecrire la matrice d’itération de la méthode itérative de Jacobi. Pour quelles valeurs de 
a cette méthode converge-t-elle ? 
-Ecrire la matrice d’itération de la méthode itérative de Gauss-Seidel. Pour quelles valeurs de 
a cette méthode converge-t-elle ? 


3 0 =T 
A=| 2 4 2 
—1 0 3 


Déterminer le polynôme caractéristique de A. 

Démontrer que A est diagonalisable et déterminer une matrice D diagonale et une matrice P 
inversible telles A = PDP) 

Calculer A"! . A3 parla méthode de Hamilton avec résultats numériques. 


Exercice 3 (Tpts) 


Soit A la matrice suivante 
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Examen de Rattrapage 
Méthode Numériques 
2016 


Questions de cours (8pts) 


a) Démontrer l’interpolation polynomiale de Lagrange d’ordre n donné par la formule 
suivante : 


n 
Pn(r) = > ) Lit) fi 
i—=0O 


n 


Ln (T) e II 


j=0. j#i Vi 


LT — Tij 


— Tj 


Avec 


b) Le choix d’une des méthodes de résolution des systèmes non linéaires est 
conditionné par quels critères ? 


Exercice 1 (6pts) 


- Résoudre par la méthode de Gauss Jordon en détail le système linéaire suivant et 


déduire le déterminant et la matrice inverse 


x+z=2 
—-x+y+t=1 
—X—y+z+t=0 
—x—y-zZz+t= -2 
Exercice 2 (6pts) 


1/2 —1/2 3/2 
On considère la matrice suivante A= |—1/2 1/2 —1/2 
0 0 1 
Coché la ou les bonnes réponses 
La matrice A a deux valeurs propres distinctes. | | 
La matrice À a trois valeurs propres distinctes. | | 


La somme des valeurs propres de A est égale à leur produit. | | 


La matrice À est diagonalisable. | | 
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Examen Final 
Méthode Numériques 
2017 
Exercice 1 (6pts): 
Soit le système linéaire suivant : 


2x1 + X2 + 4X4 = 2 
— 4x; = 2x2 + 3x3 Ka? 7x4 = —9 
4x1 + X2 — 2X3 + 8x4 = 2 
EK — 12x3 — X4 = 2 


-Résoudre le système AX = b par la méthode de Gauss 

-Déduire le déterminant 

-Extraire la deuxième ligne, la troisième colonne el la diagonale de la matrice A par commande 
Matlab avec résultats. 

-Donner les résultats des instructions Matlab : D= diag (A) ; C= A(:,1) ; R=([C_DY -A(1 :2, :)) 


Exercice 2 (6pts): 


Dans l’espace vectoriel R? , on considère l’endomorphisme U dont la matrice dans la base 


canonique est la matrice A 
3 L-=1 
A=|0 2 0 
1 1 1 


-Calculer le polynôme caractéristique de la matrice A. 

- Calculer les valeurs propres 

-En déduire le déterminant etla trace par le théorème des valeurs propres 

-Calculer la solution détaillé de A? et AT (par le Théorème de Caylay Hamilton avec 
résultats numériques) 


Exercice 3 (8pts): 

Soit la fonction suivante f(x) = x — cos (=) = 0 

-Séparer graphiquement les racines de la fonction dans l'intervalle [0, 4]. 

-Montrer que la fonction f(x) admet une racine unique a sur l'intervalle [0, 1] 

-Calculer les quatre premières itérations avec quatre chiffres significatifs après la virgule par la 
recherche dichotomique sur l'intervalle [0,1] 

-Écrire la suite de Newton associée à f(x) dans l'intervalle [0, 1] avec démonstration. 

- Quelle sont les conditions d'application de la méthode de Newton. Vérifier si elles sont 
applicables. 

- Calculer les quatre premières itérations par Newton avec zu = 1. Conclure. 
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Examen de Rattrapage 
Méthode Numériques 
2017 
Exercice 1 (6pts): 
Soit le système linéaire suivant : 


2x+y—4z=8 
3x + 3y — 5z = 14 
4x + 5y — 2z = 16 


-Résoudre le système AX = b par la méthode de Gauss Jordon ainsi que sa matrice inverse 471 
-Déduire le déterminant en utilisant la matrice triangulaire 
-Résoudre le même système en utilisant la méthode de cramer ? conclure ? 


Exercice 2 (6pts): 


-Donner avec démonstration la formule matricielle de Jacobi a partir du point fixe 

- Démontrer la suite approchée de la sécante pour résoudre l'équation non linéaire de type 
f @) = 0 

- Du point de vue résolution et convergence, qu’elle est la meilleur méthode Gauss, Jacobi, 
Gauss Seidel ou Gauss Jordon ? pour quoi ? 

- Dans la résolution des systèmes linéaires, il existe deux types de méthodes: directes et 
itératives, quelle est la meilleur méthode ? pour quoi ? 


Exercice 3 (8pts): 
1) Résoudre par la méthode de Jacobi le système suivant 


10x + y +z = 12 
2x + 10y +z = 13 
2x + 2y + 107 = 14 


On prendra pour approximation initiale x® = (1.2, 0 , 0) et on cherchera une approximation de 
la solution avec une précision de 1072 


2) Calculer le rayon spectral de la matrice d'itération de Jacobi ? Conclure 
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Examen Final 
Méthode Numériques 
Mai 2018 


Exercice 1 (7 pts) 


Soit a un paramètre réel et soient les matrices A et le vecteur b définis par 


5 2 —1 6 
A=|2 6 -3 b=|4 
a 1 4 7 


À quelle condition sur a, la matrice A est inversible (admet une matrice inverse) ? 
Soit a=- 1. Résoudre le système linéaire AX = b par la méthode du pivot de Gauss 
Déduire le déterminant 

Calculer les valeurs propres de À 

Déduire la trace et vérifier le déterminant de A 


SE RS MEL E 


Exercice 2 (6 pts) 
Soit f une application de R dans R définie par 


f(x) ei — 4x2 
1. Montrer qu’il y a une racine a comprise dans l’intervalle [0 1]. 
2. Trouver les approximations de la racine (2 chiffres après la virgule) en utilisant la 
méthode de Newton pour la recherche des zéros de la fonction f avec nz ` en 3 


itérations. 
3. Discuter sur la convergence des méthodes de résolutions des systèmes d’équations non 
linéaires. 
Exercice 3 Cours _(7 pts) 


Considérons le système linéaire AX = b avec 


a 0 y 
E a ` 
0 ô a 


Avec a, p, ô et y sont des paramètres réels de la matrice A. 
1. Donner des conditions suffisantes sur les coefficients de A pour avoir une convergence 
de la méthode de Jacobi et la méthode de Gauss-Seidel. 
2. Ecrire les formules de Jacobi etde Gauss —Seidel au avec démonstration matricielle. 
3. Donner les matrices D, U et L qui correspondent à la décomposition A = D — U — L de la 
méthode de Jacobi 
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Implémentation sous Matlab de la méthode de gauss pour un système d'ordre 3 


clc, clear all, close all 

%Nombre de variables et d'équations 
n=3; 

%lnitialisation de la matrice du système 
a=[4 -22;54-3 -2;2 3 -1|; 
%lnitialisation du vecteur b 

b=[-2 -11 -1]'; 

%Formation de la matrice augmentée 
A=[a b]; 


%Algorithme de triangularisation de Gauss 
for k=1:n-1 
for i=k+1:n 
A(i,:)=A(i,:)-(A(i,k)/A(k,k))*A(k,:); 
end 
end 
A 
%Extraction de la solution du système d'équations 
for i=n:-1:1 
s=0; 
for j=i+1:n 
s=s+A(i,j)*x(); 
end 
x(i)=(A(i,n+1)-s)/A(i,i); 
end 
Kë 


2 3 -1 -1 


0O 0 -18 -36 
Après substitution on trouve la solution suivante 


X=(-1 ,1,2) 
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